Discrete time quantum walk with nitrogen- vacancy centers 
in diamond coupled to a superconducting flux qubit 
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We propose a quantum-electrodynamics scheme for implementing the discrete-time, coined quan- 
tum walk with the walker corresponding to the phase degree of freedom for a quasi-magnon field 
realized in an ensemble of nitrogen- vacancy centres in diamond. The coin is realized as a supercon- 
ducting flux qubit. Our scheme improves on an existing proposal for implementing quantum walks 
in cavity quantum electrodynamics by removing the cumbersome requirement of varying drive-pulse 
durations according to mean quasiparticle number. Our improvement is relevant to all indirect-coin- 
flip cavity quantum-electrodynamics realizations of quantum walks. Our numerical analysis shows 
that this scheme can realize a discrete quantum walk under realistic conditions. 

PACS numbers: 05.40.Fb, 03.67.Lx, 03.65.Yz 



I. INTRODUCTION 

Spin ensembles that are coupled to quantum electrody- 
namics (QED) systems and solid-state qubits are promis- 
ing hybrid architectures [l| for quantum memories [2^ 
and for distributed quantum networking [3]. Nitrogen- 
vacancy (NV) centers in diamond are outstanding among 
the solid-state spin ensembles Q due to their long spin 
relaxation and dephasing times, scalability [5[ and ex- 
istence of well-established techniques for their manipu- 
lation NV center quantum memory has recently 
been demonstrated Q. 

Quantum memories require fast and reliable quantum 
state transfer, which is closely related to the quantum 
walk (QW) problem [1QI - [12|. I n particular, the discrete 
time QW (DTQW) [^^[ilKlsl has been highlighted for 
universal quantum computation [16] and for quantum 
state transfer [17]. NV centers that are coupled to super- 
conducting transmission line resonators are considered 
before for continuous-time QW [l8|. Here we propose a 
scheme for implementing DTQW with NV centers. 

Our scheme is based upon the proposals for DTQW 
in phase space introduced for cavity QED (CQED) fiol - 
[21] and circuit QED (cirQED) [22]. Like some other in- 
triguing proposals [23-25], QWs in CQED and cirQED 
schemes have not been demonstrated yet; despite some 
successful realizations in different systems [26M3l(. Our 
procedure could lead to the first solid-state realization of 
the proposed DTQW in phase space. 

We consider an NV center ensemble coupled to a su- 
perconducting flux qubit [32] , which plays the role of the 
coin in DTQW. Weak excitations of NV centers can be 
described in terms of bosonic quasiparticles. We take into 
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account all the crystallographic classes of the NV centers 
so that the single mode coupled to the qubit becomes an 
effective privileged mode, which is a particular superpo- 
sition of all the modes corresponding to the quasiparticle 
excitations in four different NV center crystallographic 
classes. We derive an exact bilinear spin-boson model 
that reduces to Jaynes-Cummings model (JCM), same 
with that in CQED or cirQED settings, under paramet- 
ric approximation. 

We employ the similar but improved DTQW scheme 
proposed for JCM to our case, by removing the cumber- 
some requirement of varying drive-pulse durations. This 
allows for the direct realization of DTQW with regu- 
lar time steps and simplifies the implementation in all 
indirect-coin-flip cavity QED realizations of QW includ- 
ing our full solid-state proposal. 

This work is organized as follows. Sec. |TT] presents the 
model Hamiltonian, describing the coupling of the NV 
center excitations to the flux qubit. The QW Hamilto- 
nian is obtained by employing canonical transformations. 
Sec, mil describes how the model Hamiltonian is used to 
evolve the initial state of the QW using the quantum 
master equation by taking into account the decay chan- 
nels from the dephasing and relaxation of the qubit as 
well as the relaxation of the NV centers. The dynamics 
of the mean number of NV center quasiparticles (called 
quasimagnons), the populations of the qubit levels rep- 
resenting the coin operation, and the Holevo standard 
deviation are reported. The phase distribution and the 
Wigner function plots are also provided to reveal details 
of the QW dynamics in the quasimagnon phase space. 
A linear fit analysis for the Holevo measure is presented 
and used to distinguish a QW from the classical random 
walk evolution. Sec. |lVl is devoted to the summary and 
conclusions. 
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II. MODEL SYSTEM 



B. Flux qubit Hamiltonian 



A. NV center ensemble Hamiltonian 

A theoretical model of NV centers in diamond can be 
developed based on the six-electron description of the 
electron spin triplet ground state of NV~ jssi, . The 
hyperfine ~ 2 MHz and quadrupole ~ 5 MHz interactions 
are negligible [35] relative to zero-field splitting ~ 2.88 
GHz. The Hamiltonian of a single NV center can be 
written as 



(1) 



where S is the spin-1 matrix. The zero-field splitting 
matrix D becomes diagonal in the molecular body frame; 
as a result, the NV center Hamiltonian becomes 



H = D 



(2) 



where S = 1 and D /2it = 2.87 GHz [sg] are now scalar. 
By taking the quantization axis as the alignment of the 
NV center, the Hamiltonian can be expressed in the spin 
basis \S^m) = \m) as 



H = D 



\m){m\ 



(3) 



There are four possible quantization axes in the NV 
center, which are all equivalent. We need to take them 
into account when we consider an ensemble of NV cen- 
ters. Let us denote these axes by / = 1, . . . , 4 and intro- 
duce collective operators, which are analogs of Hubbard 
operators for strongly correlated electron systems [s^, 
such that 



Nj 



(/) 



(4) 



Here, | m)j is the spin state of the j-th NV center aligned 
along the crystallographic (molecular) axis-/ which is 
taken as the quantization axis. Nf denote the number of 
NV centers aligned along the quantization axis /. Oper- 
ators XJ^^ can be conveniently related to bosonic oper- 
ators Xfm such that 



in terms of which the Hamiltonian becomes 

mf 



(5) 



(6) 



The operators Xfm play the role of Schwinger bosons of 
SU{Nf) subject to constaint 



Nf |38hl^ 



A fiux qubit consists of three or more superconduct- 
ing Josephson junctions, one of which is smaller than the 
others [4l|, |42| . The phase quantization condition for the 
superconducting wave functions around the qubit loop 
reduces the system to a two-dimensional phase space, 
where the qubit can be effectively considered as a two- 
particle system subject to a double- well potential 0. A 
two-level approximation can be made to the band struc- 
ture of the potential by assuming that the external fiux is 
in the vicinity of half of the fiux quantum, where the two 
lowest bands are closest to each other. A small overlap 
of the wave functions of the two levels opens the qubit 
gap Aq, which can be tuned [44l446| by replacing the 
smallest Josephson junction with two parallel junctions, 
a DC-SQUID. On the basis of the two counter-rotating 
persistent currents in the qubit loop, the qubit Hamilto- 
nian can be written as 



q - 2 i^q^^ 



(7) 



where Cq is the magnetic energy bias of the qubit that 
contributes away from the symmetry point, when there 
is a net persistent current in the qubit loop. 

Pauli spin operators (Jx,y,z can be conveniently trans- 
formed from the bare current basis to the dressed current 
basis by U = exp {—iOay) with = arctan (Ag/e^). The 
Hamiltonian becomes 



^q = 



(8) 



with the qubit frequency uOq = J (A^ + e^)/2. 



C. Flux qubit-NV ensemble coupling 

We assume the same conditions as in a recent experi- 
ment [13 , wherein the fiux qubit was coupled to the NV 
center ensemble located at the top of the central region 
of the fiux qubit loop. The dimensions of the ensem- 
ble are sufficiently small for the spatial variations in the 
magnetic field over the ensemble to be ignored. In such 
a case, the magnetic fields generated by the super cur- 
rents are expected to interact with the ensemble homo- 
geneously. In the persistent current basis the electronic 
Zeeman interaction is given by 



E 



(9) 



where are the total spin-1 operators for the NV en- 
semble that are related to the generators of the SU(3). 
Nuclear Zeeman interaction is three orders of magnitude 
smaller and neglected. 

The coupling strengths ^£ = 2^eMB^/a5 where Qe is 
the electronic Lande g factor, /iB is the Bohr magneton 
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and Bfa is the a-component of the magnetic field acting 
in the molecular frame of the crystallographic class /, can 
be adjusted by the orientation of the diamond sample. In 
a typical experimental setting, the crystal axis (001) is 
parallel to the magnetic field, which makes = the 
same for all crystallographic classes / and gz becomes 
negligible j32|, [Z^- In addition, the y direction can be 
selected to make gy relatively smaller than g^ = g [47]. 
These lead to the reduced interaction of the form 

i^int = ^^^.(X^VX+Vh.c.), (10) 

/ 

= ^^cr^(xJ._x/o + 2;J.^^/o + H.C.), (11) 

/ 

where H.c. stands for Hermitian conjugate. 

We shall proceed under the assumptions of weak ex- 
citations and a large number of spins in the ensemble. 
Experiments can achieve coupling of N ~ 3.1 x 10^ NV 
centers directly above the fiux qubit [13 • This number is 
sufficiently large for our weak-excitation conditions. We 
consider about 5 — 10 excitations in our simulations. The 
quasiparticles of NV center excitations are analogous to 
quantum particles associated with the spin waves, and 
hence, we call them quasimagnons. We introduce quasi- 
magnon operators as follows: 

X/+ ^ a/, Xf- bf, (12) 

and proceed by replacing Xfo ^ \/^f which is analo- 
gous to the parametric pump approximation in quantum 
optics Parametric approximation has quite a wide 
range of validity, which is especially enhanced for initially 
coherent states [49]. To the first order the error is as 
large as the of the square root of the mean quasimagnon 
number. Fractional error, relative to the mean number 
of quasimagnons, is about the size of relative fiuctuation 
which is negligible in our case \fN /A/" ~ 0.01. 

We consider initially coherent spin ensembles for our 
quantum walk procedure as well, which could make our 
procedure potentially applicable for higher number of ex- 
citations, or equivalent ly for smaller number of NV cen- 
ters coupled to the fiux qubit. A systematic study of how 
tolerant our procedure to smaller number of NV centers 
however needs further examinations without parametric 
approximation, which is beyond the scope of our present 
investigations. 

The Hamiltonian of the NV centers then reduces to 

i^NV = D{a\af + b\bf). (13) 
/ 

We interpret the presence of two quasimagnon operators 
by associating them with quasispin waves with two possi- 
ble polarizations. The final Hamiltonian is then an eight- 
mode model, in which there are four quasimagnon modes, 
each with two different polarization degrees of freedom. 

By employing the rotating wave approximation for 
g <C D^uOq^ transforming to the dressed current basis. 



and using the quasimagnon picture, we can express the 
interaction as 

^int = ^[^/(^+(a/ + hf) + a_(4 + b])l (14) 
/ 

with Tjf = g\/Nf being the collective coupling strengths 
for the quasimagnons and ffux qubit system. 

This Hamiltonian can be simplified by distinguishing 
sets of "coupled modes" and "uncoupled modes," anal- 
ogous to Morris-Shore transformations [50], which are 
combinations of quasimagnon modes such that they are 
respectively interacting and noninteracting with the fiux 
qubit. We list these modes in the Appendix. Eventually, 
only a single effective privileged mode becomes coupled 
to the qubit. The annihilation operator of that mode is 
given by 

^ = 7f^Ev^(«/ + M- (15) 

The interaction Hamiltonian is then simplified to the 
usual Jaynes-Cummings model, though the single priv- 
ileged (bright) mode is actually a spesific superposition 
of all the quasimagnon modes. 

Hint = r]{(Jj^c + (j-c^), (16) 

where the collective spin-qubit coupling strength is = 
V2Ng. 

D. Quantum walk Hamiltonian 

We introduce an additional term that represents the 
driving field on the diamond sample 

Hd = e{t){ce'^^' + c^e-^^^'), (17) 

where ujd is the drive frequency and e{t) is the time- 
dependent amplitude of the drive pulse. We assume that 
the drive is of square- wave form with pulse duration tn 
and amplitude eo- The total Hamiltonian of the qubit 
ensemble becomes the Jaynes-Cummings model Hamil- 
tonian such that 

H = '^az^Dc^c^r]{a^c^a-C^)^Hd{t). (18) 

This Hamiltonian is the same as that proposed for a 
QW in cirQED, although in our case, the role played 
by the microwave photons is assumed by the quasi- 
magnons. Here, the qubit-NV center ensemble coupling 
T] is about r]/27r ^ 70 MHz ^]. In contrast to sponta- 
neous decay and the dephasing of the qubit, the qubit- 
NV center ensemble interaction is in the strong-coupling 
regime. By employing a rotating frame transformation 
Hr WH'U - K with U exp {-iAt), where we take 
A = ujd{c^c -\- the Hamiltonian can be written in 

the form 

H = -Scc'^c - SdY + ^(c^^- + ca+) + e(c + c^). (19) 
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Here, 5d = ujd — and 5c = ujd — D are introduced. 

The detuning between the qubit and the NV center 
ensemble is A = dd — Sc = D — Uq, which can be adjusted 
by ujq to satisfy A ^ 77 so that the Frochhch canoni- 
cal transformation [51-53] H' = exp {—S)H exjp (S)^ with 
the transformation operator 



S 



V 



■Cr+C), 



(20) 



can be applied to obtain an effective Hamiltonian 



X 



SrC^C - 



(21) 



where I2 is the 2x2 unit matrix; and the parameters 



27?e(t) 



(22) 



are associated with the conditional rotational shift and 
the Hadamard coin terms in the Hamiltonian. The 
Hadamard coin term only operates when the pulse is on, 
whereas the shift operation is constantly active. 

We note that in the QW proposal for cirQED, the pa- 
rameter Qji is found to depend on uj^ by the relation 
Qr = 2r]e{t)/5c. The dependence of Qr on the drive 
frequency results in a QW with varying step sizes in 
time, and complicated pre- adjustments in pulse duration 
are therefore required. Here, we do not have a drive- 
frequency-dependent coin operation and we have a fixed 
step size as well as fixed pulse duration. This makes the 
implementation of the QW more straightforward. In con- 
trast to our approach described here, a displaced frame 
transformation relative to the resonator field was previ- 
ously used in the cirQED model of a bit-fiip gate [H^ . 
in addition to the Frochlich transformation. This lead to 
vanishingly small cavity photon numbers and negligible 
ac Stark shifts. In our approach, we do not use a dis- 
placed frame, which allows for the presence of a larger 
number of quasimagnons. 



III. RESULTS AND DISCUSSION 

The dynamics of the effective Hamiltonian is subject 
to decoherence processes. By denoting the qubit dephas- 
ing rate as 7^, the qubit relaxation rate as 71, and the 
decay of the NV center excitations as F, we can write 
the quantum master equation of the system, under the 
Markovian and Born-Markov approximations, which are 
suitable for our weakly coupled system as 

p = -i[H, p] + TV[c]p + jiT)[a.]p + ^V[a,]p. (23) 

Decoherence processes associated with an operator x = 
c, cr_ , az are described by the Liouvillian superoperators 
V[x]p = {2xpx^ — x^xp — px^x)/2 in Lindblad form. We 



introduced the effective Hamiltonian to illustrate the un- 
derlying QW dynamics, but obtain quantitative results 
from the Hamiltonian of the system given in Eq. (p!9|) . 

The effective Hamiltonian is in We do not simulate the 
effective Hamiltonian to avoid further approximations. 
The effective Hamiltonian only serves to interpret and 
design the dynamics of the system. 

We^olve the master equation using the QuTip pack- 
age [56| in PYTHON software. The initial condition is 
considered to be 



\m) = Hi\9)+i\e))/V2, 



(24) 



where the quasimagnons are in a coherent state with am- 
plitude a and the qubit is in a superposition of its ground 
and excited states. It is assumed that excitations of the 
NV center ensemble is initially prepared in a coherent 
state. We consider a Fock space of quasimagnons with 
the dimension 17, which is sufficiently large to accom- 
modate the selected a = 3 in the numerical analysis. 
A classical magnetic field pulse can be used to prepare 
the quasimagnon coherent state analogous to the typical 
driven oscillator case [57], as for large number of NV cen- 
ters the atomic coherent states converges to the bosonic 
coherent states [58]. The modes to be included in the 
initial coherent state could be selected by external bias 
magnetic fields. The initial coherent superposition of the 
fiux qubit levels can be generated by a resonant magnetic 
pulse of area tt. 

Most of the parameters we use in our simulations can 
be taken as the same values with those used in the orig- 
inal circuit QED and cavity QED proposals [21, 22] ex- 
cept 7i,r,I) and eo- Flux qubit has a wide range of 
tunability and we set its gap frequency as ujq /27r = 7 
GHz [42|- Recently developed fiux qubits can have de- 
phasing times in the order of few microseconds [59| and 
we set 70/27r = 0.31 MHz. Recent experiments allow for 
direct coupling between the fiux qubits and the NV cen- 
ters with strengths of about ~ 70 MHz [47^. Number of 
NV centers and the fiux qubit size can be used to con- 
trol the strength of the coupling, which is taken here as 
r]/27r = 100 MHz. 

As our improved procedure is general and applicable 
to previous circuit QED or cavity QED settings [2l|, |22| , 
we shall first use the original values 7i/27r = 0.02 MHz 
and eo /27r = 1 GHz for comparison of our improved pro- 
cedure with the original ones in Refs. [2l|,[22|. The minor 
difference of no consequence is that we use I)/27r = 2.87 
GHz [sil for the zero-field splitting between the | m = 0) 
and I m = ±1) ground state spin triplet levels of the 
NV center. This corresponds to cavity or transmission 
line mode frequency, which is taken as 5 MHz in Refs. 
[2l|,[23. 

Following the comparison, we present the results for 
the actual fiux qubit and NV ensemble parameters. Re- 
laxation rate of the fiux qubit is of the same order with 
its dephasing rate so that we set 7i/27r = 0.31 MHz [59|. 
Longitudinal relaxation time Ti of the NV center elec- 
tron spin is in the order of few milliseconds at room tem- 
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FIG. 1: (Color online) (a) Time (t) dependence of the mean number of quasimagnons {uc) and (b) of the qubit-level populations 
Pe,g for the initial state |'0(O)) — \a — 3)(|^)+i|e))/v^- (^c) fluctuates in time, mainly when the pulse is on (coin toss). Regular 
intervals of the pulse in (b) indicate the regular discrete time steps of the QW. The parameters used in our simulations are 
cjg/27r = 7 GHz, 0/271 = 2.87 GHz, 7y/27r = 100 MHz, eo/27r = 1 GHz, 7i/27r = 0.02 MHz, 70/27r = 0.31 MHz, and V/2ti = 0.1 
MHz. The repetition rate of the driving pulsed laser is taken as ujp = with d = 16 being the number of sites on the cycle 
of a QW. 



perature [s^ and gets longer by reducing the tempera- 
ture. The dephasing time of the spin T2 is in the order of 
0.3 — 2 ms and gets longer by increasing the purity of the 
diamond sample isotropically [60]. In the case of quasi- 
magnons, the decoherence rate F would be collectively 
enhanced, similar to the collectively enhanced coupling 
strength. For = 3.1 x 10^ and T2 = 2 ms, we find 
F/27r = 2.78 which is larger than the circuit QED and 
cavity QED cases where F/27r = 0.1 MHz [U, [22|. To 
compensate it we take larger drive strength eo/27r = 10 
GHz. 

Relative to the decoherence ^ and dephasing [62] of 
the qubit, as well as to the relaxation of the NV cen- 
ters, the system is in the strong-coupling regime, though 
still not in the ultrastrong-coupling regime, which would 
challenge the validity of rotating wave and Born-Markov 
approximations. 

The repetition rate of the driving pulsed laser is taken 
to be ujp = with d = 16 being the number of sites on 
the cycle of a QW. This yields the duration of each time 
step of the QW as tp = 27: /uop. This time covers the pulse 
duration of tu = 7r/(A/2^i?) as well as the time between 
pulses. Hadamard coin operation induced by the pulse 
can be described by the Hamiltonian 



He, 



(25) 



if we take the operation time as tn = 7r/{V2QR); and if 
we decouple the quasi-magnon dynamics from the qubit 
dynamics. For that aim, let us drop the constant term 
(x/2)l2 and rewrite the effective QW Hamiltonian in 
Eq. (jni) as 



^ {xc^c-n)a^ - ScC^c^e{c^c^), (26) 



under the condition of 



Xn 



X 
2 



2 



(27) 



This condition can be satisfied by taking the pulse fre- 
quency to be 



Ud = (2n + l)x + cjg - 



(28) 



^eff = ^c< 



where n = |ap is the initial number of quasi-magnons. eo 
is the amplitude of the square pulse. If the number fluc- 
tuations are sufficiently small, then the dynamics of the 
quasi-magnons are uncoupled from the qubit degrees of 
freedom, whose evolution would be a faithful representa- 
tion of that of an ideal Hadamard coin. Quantitatively, if 
we denote the number fluctuations by (5n, the condition of 
validity can be expressed as x^^ <^ ^i?/2 or Sn <C eo/r], 
which becomes (5n <C 10 — 100 for our parameters. Our 
simulations indicate that number fluctuations remain in 
the order of few quasi-magnons; and thus the coin oper- 
ation can be simulated. In contrast to the cirQED and 
CQED QW proposals, we do not change ujd by changing 
n at every step. Furthermore, the pulse duration, and 
hence, the time steps, are fixed in our approach. These 
two differences allow for simpler implementation as well 
as for direct interpretation in terms of DTQW with reg- 
ular time steps. 

We first examine the behavior of the number of quasi- 
magnons, {ric) = (c^c), and the qubit excited and ground 
level populations, which are respectively defined to be 
Pe = (l + (cr^))/2, andP^ = (l-(cr^))/2. Their time evo- 
lutions are shown in Fig. [H The initial number of quasi- 
magnons is nc(0) = 9. Fig. [Tal indicates that during the 
course of dynamics mean number of quasimagnons fluc- 
tuates because of the drive and dissipation. Relatively 
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Strong fluctuations in the mean quasimagnon number oc- 
cur mainly when the pulse is on, or when the Hadamard 
coin is tossed. When the pulse is off, a conditional phase 
shift operation characterized by the parameter x is set 
into action. During this rotational QW step in quasi- 
magnon phase space, fluctuations of the (ric) are much 
smaller. The dynamics of the coin, or the populations of 
the qubit levels, can be followed by the Fig. flbl We have 
marked the regions during which the pulse is on. The dy- 
namics and exchanges between the upper and lower level 
populations reveal that the qubit evolves in close cor- 
relation with the simulation of a pure Hadamard coin. 
Exchanges only occur when the pulse is on, or at the end 
of the coin toss. 




0.1 0.2 0.3 

t(^s) 

FIG. 2: (Color online) Time (t) dependence of the Holevo 
standard deviation measure an in the system for the initial 
state \ip(0)) — \a = 3) {\g) -\- i\e)) / ^2. Its decrease beyond the 
seventh step is due to the cyclic wrapping-up efl'ect in terms 
of the definition of the Holevo measure; it does not represent 
the actual behavior of the QW. The parameters are the same 
as Fig. □ 

QW in the quasimagnon phase space needs to be char- 
acterized by examination of the dispersion of the angu- 
lar variable (j). For cyclic variables the variance is con- 
veniently determined by the Holevo standard deviation 
which is deflned as [63| 

<TH = V^l (e'^) -1, (29) 

where 

I (e'*) 1 = 1 r d0P(0)ei* I (30) 

J —TV 

is the sharpness of the phase distribution. Sharpness is 
bounded between 0, for spiked distribution at the mean 
phase, and 1, for a flat distribution. It can be interpreted 
as the success or fldelity of the phase distribution to es- 
timate the mean phase [63, Ell . Here we take the mean 
phase as and choose the domain as [— 7r,7r]. The phase 
distribution P(0) is determined from the reduced density 
matrix of the quasimagnons pm by using 

P{<P) = lim {4>\pm\ 4>), (31) 



2 




log(t) 



FIG. 3: (Color online) A linear fit to the Holevo measure cth 
in Fig. [21 A logarithmic scale is used for both axes. The 
initial starting point in time is taken as zero and not used on 
the logarithmic scale. The slope of the line is calculated as 
^^0.96, which is in the ballistic regime and higher than the 
corresponding value of 0.5 for a classical random walk. The 
parameters are the same as Fig. [T] except that the time t here 
is dimensionless, scaled by 1 /xs. 



where the phase states are deflned to be 

M-l 

\'^)=mY. e'""* I (32) 

M=0 

The quasimagnon reduced density matrix is calculated 
by, Pm = Trq(p), taking the trace over the qubit de- 
grees of freedom of the density matrix of the system p, 
which is found by solving Eq. (|23|) . If the phase vari- 
ance (A0)^ = (0^) - (^)^ is smah ((A^)^ <C 1) then an 
becomes equivalent to the A0. In contrast to localized 
distributions, Holevo variance is inflnitely large for flat 
distributions with zero sharpness. In our case, this sit- 
uation corresponds to large phase diffusion due to deco- 
herence in the system. We examine the Wigner function 
evolution to determine the number of steps for which the 
Holevo measure can be used for reliable information on 
the phase variance. 

Fig. [2] shows that (Jh increases in an almost smooth 
manner during the conditional phase shifts, whereas dur- 
ing the coin toss operation, it shows rapid oscillations. At 
the eighth step, more than half of the cycle is covered by 
the Wigner function (see Fig. [5]). The Holevo measure 
becomes unsuitable for describing the second order phase 
variance beyond this point. Similar behavior is encoun- 
tered in the cirQED QW proposal [21]. To check whether 
the phase variance is faster than that of the classical ran- 
dom walk, we employ a linear flt analysis for the flrst 7 
steps using MATLAB. The result is shown in Fig.[3l The 
upper panel in Fig. [3] reveals that there is a linear rela- 
tion between the number of regular time steps and the 
logarithm of the Holevo standard deviation. The slope is 
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determined to be ~0.96 ±0.11, where 0.11 is the hnear- 
regression standard deviation. This result is close enough 
to the ballistic case to declare it to be ballistic given that 
phase diffusion due to decoherence and wrapping due to 
periodicity of phase have occured. In the case of the 
classical random walk, the slope would be 0.5. 

The observation of further QW signatures is to be ex- 
pected, similar to earlier studies on the same Hamilto- 
nian for similar parameters. For example, the phase dis- 
tributions at the end of the first four steps are reported 
in Ref. [21]. Interpretation of the phase distribution is 
not trivial because of the effects of the pump field. It is 
found that as the pump intensity grows, the distribution 
becomes asymmetrically skewed in one direction. More- 
over, some peaks become narrower. We present the phase 
distribution in Fig. SI We used M = 256 phase states in 
our numerical calculations. Phase distribution is skewed 
counterclockwise because of the pumping field, consistent 
with the observations reported in Ref. [21]. Because the 
walk occurs on a circle, the left and the right boundaries 
(j) = ±7r in Fig. |4]are the same. The initial phase distri- 
bution of the QW starts at an angular location of (j) = 0. 
The spread of the distribution is quadratically enhanced, 
relative to that of a classical random walk. 

Wigner functions at the end of the first eight steps 
are shown in Fig. [5l The spread along the angular di- 
rection as the QW progress is accompanied by the in- 
creasing deformation along the radial direction due to 
the number fluctuations. QW splits into other circles 
at different mean number of quasimagnons and dephas- 
ing occurs. Holevo measure increases beyond 1 as the 
distribution gets more and more flat due to significant 
contribution of the phase diffusion after 7 steps. Thus 
we did not use the Holevo measure to characterize the 
QW for more than 7 steps. 

We now examine the effect of decoherence rate of 
presently available NV center ensembles with r/27r = 
2.78 MHz [60] which is about 10 times larger than the 
photons in cirQED and CQED QW proposals. To com- 
pensate faster loss of quasimagnons, we take 10 times 
stronger drive with eo/27r = 10 GHz in the simulations, 
which is found to be sufficient enough to stabilize the cir- 
cular trajectory of quantum walkers in the phase space 
by maintaining (ric) ~ 6 as can be seen in Fig. [6l For 
this drive strength, the drive duration becomes 10 times 
shorter, due to the relation tn ~ 1/^i? ~ l/eo- The 
spikes in the figure correspond to the times when the 
drive is on or the action of coin toss. The Holevo stan- 
dard deviation shown in the Fig. [71 In this case phase 
diffusion becomes significant after four steps. When we 
determine the linear fit in logarithmic scale we find the 
similar behavior as in Fig. [3] with the slope 0.98 ± 0.42, 
where 0.42 is the linear-regression standard deviation for 
the first four steps. Such a few number of steps yields 
large linear-regression standard deviation error in the 
slope. The slope is nevertheless close enough to be de- 
clared to be ballistic, given that phase diffusion due to 
significantly large decoherence and wrapping due to pe- 
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FIG. 4: (Golor online) The phase distribution P = P($) 
as a function of the angular variable $ of the quasimagnon 
phase space is shown at the end of the (a) first, (b) second, (c) 
third, and the (d) fourth steps of the QW. The initial phase 
distribution, wherein the QW starts at $ = 0, is not shown. 
The parameters are the same as Fig. [Tl $ is in radians. 

riodicity of phase have occured. 

IV. CONCLUSIONS 

We proposed a fully solid-state implementation of a 
DTQW of NV centers in diamond coupled to a flux qubit. 
The QW takes place in the phase space of the quasi- 
magnon field by indirectly flipping the flux-qubit coin. 
The proposed scheme is an improved version of the one 
originally developed for cirQED and cavity QED systems. 
In our method, the time steps are of fixed duration and 
no alteration of the pulse is required. This improved 
method is applicable to general cirQED and cavity QED 
QW scenarios. We did not find any significant reduction 
in the quality of the QW, characterized by the spread of 
the Holevo standard deviation. 

If the NV center ensemble is sufficiently large and 
weakly excited, its excitations can be treated as bosonic 
quasiparticles, called quasimagnons. A driven Jaynes- 
Cummings model for the quasimagnons and the flux 
qubit system is established. Further canonical transfor- 
mation of the Jaynes-Gummings Hamiltonian to a QW 
on a cycle Hamiltonian is then applied, following a pro- 
cedure similar to that used for a QW in cirQED. Cal- 
culation of the Holevo measure, by taking into account 
the dissipation channels, and then fitting a curve to it on 
a logarithmic scale, reveal that the spread of the QW is 
quadratically faster than that of a classical random walk. 
Further, QW signatures are explored in the phase distri- 
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FIG. 5: (Color online) The Wigner function at the end of the first eight steps of the QW. The parameters are the same as 
Fig. [T] Real and imaginary parts of the coherent state parameter in the quasimagnon phase space are denoted by x and p, 
respectively. The initial Wigner function, which is a Gaussian centered at (x = 3,p = 0), is not shown. The radius of the 
circle of the QW is determined by the square root of the mean number of the quasimagnons. Phase diffusion along the angular 
direction is accompanied by the ripples along the radial direction due to the number fluctuations as the QW progress. 




FIG. 6: (Color online) Time (t) dependence of the mean 
quasimagnon number (ric) under the same conditions as in 
Fig. [U but for the case of r/27r = 2.78 MHz and drive 
eo/27r = 10 GHz. 



FIG. 7: (Color online) Time (t) dependence of the Holevo 
standard deviation measure ctr under the same conditions as 
in Fig. Elbut for the case of r/27r = 2.78 MHz and drive 
eo/27r = 10 GHz. 



Acknowledgments 



bution and the Wigner function, and features that are 
similar to those found in the cirQED and cavity QED 
cases are confirmed. 



O. E. M. acknowledges financial support from 
TUBiTAK (Grant. No. 112T049). B. C. S. acknowl- 
edges financial support from AITF, NSERC, and CI- 
FAR. P. X. acknowledges financial support from NSFC 
11004029 and 11174052. 



9 



Appendix: COUPLED AND UNCOUPLED 
QUASIMAGNON MODES WITH A FLUX QUBIT 

Two quasimagnon modes corresponding to two qua- 
sispin waves with different polarizations are nearly de- 
generate. This allows for decoupling a combination of 
them from the flux qubit. The annihilation operators for 
the coupled and uncoupled modes are Cf = (a/ -\-bf)/y/2 
and df = {ctf — &j)/\/2, respectively. These are also 
nearly degenerate, whereby subsequent decouplings are 
possible. New annihilation operators, given by 



are deflned to decouple di2 and ^34 from the qubit. When 
the procedure is applied once more, it leads to the cou- 
pling of the qubit to a mode described by the annihilation 
operator 



c = 



y/Ni + N2C12 + VN3 + N4C34 



N 



=(E + 



/27V 



(A.2) 



C12 
C34 

^34 



y/Ni + N2 
1 

VNi + N2 ' 
1 

1 

VN3 + A^4 ' 



(VTViCi + VN2C2), 



WN2C1 - V^lC2), 



(VAr3C3 + VA^4C4), 



(V7V4C3 - ViV3C4), 



(A.l) 



with N = J] J Nf as the total number of NV centers. 
The uncoupled mode is given by 



d = 



VN3 + N4C12 - VNi + N2C34 



(A.3) 
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